On the Local Structure of Topological Charge Fluctuations in QCD 
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We consider the lattice topological charge density introduced by Hasenfratz, Laliena and Nie- 
dermayer and propose its eigenmode expansion as a tool to investigate the structure of topological 
charge fluctuations in QCD. The resulting effective density is built from local chiralities studied pre- 
viously. At every order of the expansion the density exactly sums up to the global topological charge, 
and the leading term describes the maximally smooth space-time distribution of charge relevant for 
propagating light fermions. We use this framework to demonstrate our previous suggestion that the 
bulk of topological charge in QCD does not effectively appear in the form of quantized unit lumps. 
Our conclusion implies that it is unlikely that the mixing of "would-be" zeromodes associated with 
such lumps is the prevalent microscopic mechanism for spontaneous chiral symmetry breaking in 
QCD. We also present first results quantitatively characterizing the space-time behavior of effective 
densities. For coherent fluctuations contained in spherical regions we find a continuous distribution 
of associated charges essentially ending at « 0.5. 
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Some intriguing effects in QCD, such as the large 77' 
mass and the ^-dependence, are related to vacuum fluc- 
tuations of topological charge |3| ■ Understanding the lo- 
cal structure of topological charge fluctuations is thus of 
great interest for building a detailed picture of how these 
phenomena arise in terms of fundamental degrees of free- 
dom. Another important phenomenon possibly related 
to topological charge fluctuations is spontaneous chiral 
symmetry breaking (SxSB). This is based on the propo- 
sition that the effective low-energy structure of topolog- 
ical charge fluctuations in QCD is such that in a typical 
configuration most of topological charge is concentrated 
in Nl typically non-overlapping, four-dimensional space- 
time regions Ci , each containing a sign-coherent lump of 
approximately unit topological charge (i.e. the general- 
ization of the instanton liquid picture) . // true, this would 
imply an appealing microscopic explanation for the ori- 
gin of Dirac near-zero modes (2| and hence the origin of 
S^SB j^] . The connection to fermions arises by associat- 
ing with each lump Ci a localized chiral mode which 
"would be" a zeromodc in the absence of other lumps, 
and a formation of a topological subspace of low-lying 
eigenmodes ■0* ~ J2j '^ijX' i * = 1, . . . Nl- Verifying this 
"topological mixing" scenario has long been hindered by 
inherent difficulties in interpreting the local behavior of 
topological charge density for typical configurations of 
lattice-regularized theory. We have argued 0, Q that 
meaningful information can be extracted indirectly by 
studying the local chirality of low- lying modes. The flat 
behavior of the associated AT-distribution observed in the 
initial study suggested that topological mixing was not 
the origin of the near-zeromodes, and also indicated the 
lack of extended self-dual excitations jjj. However, later 
studies revealed the double-peaked behavior which 



is qualitatively consistent with both these aspects, but 
not sufficient for their demonstration 5]. Being agnostic 
about self-duality, we observed in Ref. .5] that even with 
the double-peak structure present, the patterns in low 
eigenmodes of the overlap operator are inconsistent with 
the presence of a distinctive topological subspace, and 
hence with the topological mixing scenario. In this work, 
we propose a framework generalizing the local chirality 
method, and use it to demonstrate this conclusion di- 
rectly. The basis of our approach is the topological charge 
density g^; = \ tr^^ Dx,x ;7] associated with 75-Hermitian 
D satisfying {D,^^} = Dj^D (Ginsparg- Wilson (GW) 
fermions). Note that since tr75 = 0, one can also write 
Qx = — tr75 (C — \Dx,x) with arbitrary constant C. We 
choose C — 1 (see also Ref. Q) which guarantees that 
the eigenmode expansion of qx is properly normalized 
and satisfles the index theorem for arbitrary truncation 
(see below). In this case we have 



1 



qx 



A. 



tr75(l--7^,,,) - -^(l--)c^ (1) 



where = i^x^^b'i'x is the local chirality of the mode 
with eigenvalue A. For low- lying truncation (|A| « 0) 
this is effectively the sum of individual local chirali- 
ties. Very recently, Gattringer offered the evidence for a 
large degree of self-duality in the vicinity of peaks of an 
eigenmode-filtered action density . If the extended na- 
ture of self-dual regions is confirmed, then combined with 
our conclusions one is led to a picture of inhomogeneous, 
but essentially continuous topological charge fiuctuations 
at low energy, with a significant level of self-duality in the 
vicinity of local maxima. The origin of the double-peaked 
behavior of local chirality would then be the local attrac- 
tion of spinorial components by (anti)self-dual field as 
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argued in PJ, with topological mixing not playing a role. 
In the last part of this work, we describe the first at- 
tempt to characterize the typical amounts of topological 
charge associated with coherent fluctuations, as well as 
some interesting quantitative results that emerged. 

1. In the usual discussion of topological mixing there 
is a paradox rooted in the negativity of Euclidean topo- 
logical charge density correlator at nonzero distance, 



{qix)q{0)) < 0, |a:|>0. 



(2) 



This follows from reflection positivity and the fact that 
the operator q{x) is reflection-odd One conse- 

quence is that it is impossible that in typical configura- 
tions most of topological charge is concentrated in co- 
herent four-dimensional lumps of finite physical size (e.g. 
instantons). Indeed, if this was the case and the typical 
size of such lumps was Vc, then the average correlator 
would be positive over Vc, thus contradicting (01. This 
means that lumps Ci cannot be identified in q{x). 

These considerations do not rule out topological mix- 
ing as an effective low-energy scenario for S^SB, but 
emphasize the necessity of some short-distance filtering 
when studying these issues. This can be achieved by re- 
placing the local operator q(x) by a nonlocal one. While 
the use of nonlocal operators is generally unacceptable, 
here they will only serve as a filter to smear the singular 
contact term in (0) into a finite positive core. The dom- 
inance of coherent structures is then possible. Following 
up on our reasoning in Refs. 0,07 we propose that the 
low energy truncation of the Dirac eigenmode expansion 
for q(x) (see below) can serve as the physically motivated 
nonlocal operator needed. We stress the physical moti- 
vation because, for the problem at hand, the underlying 
issue is whether the light fermion effectively feels a col- 
lection of coherent unit lumps as it moves through the 
vacuum, or there is some other dynamics governing its 
propagation. This is to be decided by the fermion. 

The desired fermion filtering can be realized starting 
from lattice-regularized theory. Among the lattice topo- 
logical charge density operators considered, the ones as- 
sociated with GW fermionic kernels are perhaps theoret- 
ically most appealing 0. These are constrained so that 
they sum up to integer global charge defined by counting 
the exact zeromodes of the GW kernel used. The contin- 
uum gauge-fermion correspondence (index theorem) for 
smooth gauge fields is thus extended to the lattice by con- 
struction 7J. Before defining the filtered densities start- 
ing from it is useful to recall that the spectrum of 
D contains A'o zeromodes with global chirality -1-1 {Nq) 
or —1 (Nq), N2 modes at A = 2 with global chirality 
-1-1 {N2) or —1 {N2), and Np pairs of complex modes 
with zero global chirality. The dimension of Dirac space 
is = A^o + ^2 + 2 * A^p, and the topological charge is 

Q = J:.q^ = No -K = 

with q^ the set of related densities 
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;s eigenpairs 


£1 


6.00 


0.093 


14 


12 
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£2 


6.20 


0.068 


20' 
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£3 


6.55 


0.042 


32' 


5 
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£4 


5.91 


0.110 


12' 


6 
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£5 


6.20 


0.068 


20' 
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10 



TABLE I: Ensembles of Wilson gauge configurations. 
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FIG. 1: Top: Cq, C 



(1) 



and Cg(i6) 



for configuration C2 (Q = 0) 



from £a. Left: Roughness G'*' of q^''^ for the same configu- 
ration. For the full density, G'^f^ = 0.87. Right: Size of 
the positive core of the average C^(2) correlator for ensembles 
£i-£3. Sommer parameter was used to set the scale. 

representing the truncated eigenmode sum that includes 
all zeromodes and the k lowest-lying complex pairs. 
Complex eigenvalues Xj enter the sum ordered on the up- 
per (lower) half of the spectral circle. Eigenmode-filtered 



densities g^*"'' are real and q. 
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Q = - A+ thus 



leading to identical global fiuctuations satisfying the in- 
dex theorem. We have gi^'''' = qx and since q^^^ can 
be identically zero for Q — Q configurations, we consider 
g*-^-* to be the leading order in the expansion. The in- 
frared eigenmodes of D are significantly smoother than 
the underlying gauge field 0, and q^^^ is expected to 
be maximally smooth. As k increases, g'-'^-' gradually in- 
corporates more short-distance structure. This does not 
mean that there is a strict new cutoff present in filtered 
densities. However, ultraviolet fluctuations irrelevant for 
propagation of light quarks are filtered out. 

We now use the overlap Dirac operator |0 (see details 
in Ref. i5|) on Wilson gauge backgrounds (see Table P) 
to demonstrate the basic properties of filtered densities 
discussed above. We have calculated the full density q 
for configuration C2 from ensemble £4^. In Fig. ^ we show 
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its correlator Cq{r) (normalized at the origin) and com- 
pare it to C^(i) and C^(i6) for filtered densities. Note 
that Cq has a very short-ranged positive core as expected 
from IpJ. This provides indirect evidence that the local- 
ity of the non-ultralocal operator q is quite good. As for 
the correlators C^(fc) , their range is evidently larger. The 
shape stabilizes at about k = 6 and changes very slowly 
from then on. To further characterize the roughness of 
filtered densities we calculate 

a;, /J, 

and plot G^^) as a function of k. As expected, q^'^^ be- 
comes rougher as k increases. While the physical size 
of the positive core in Cq is expected to go to zero in 
the continuum limit, this is not necessarily so for Cq(k) . 
To see that, we have calculated Cq{2) for ensembles £1-^3. 
The size of the positive core in the average correlator was 
determined as < r > over the probability distribution 
given by C^(2) (r) in the range from zero up to the max- 
imal distance where (7^(2) (r) is manifestly positive (with 
errors taken into account). As can be seen in Fig. the 
size scales well, indicating that if the dominance of coher- 
ent four-dimensional structures in q^^^ can be established, 
these structures can survive the continuum limit. 

2. Introduction of fermion- filtered densities allows 
us to study the relevance of topological mixing consis- 
tently. If the topological subspace of low-lying modes 

w Q'iiX"' J i = Ij • ■ ■ I A^L is formed, then 

Nl 

xeCi j=i 

(5) 

This implies that the structure of quantized unit lumps 
Ci will be revealed in g^'^-' when all modes belonging 
to the topological subspace are included. Indeed, from 
equation |(SJ) we have 

^ gi'^) = « ±1 , 2k^NL-No (6) 

since Re A « for modes in topological subspace. More- 
over, the value of pure gauge topological susceptibility 
(« 1 fm^**) constrains Ni^ in lump-dominated configura- 
tions of volume V fm'' to be Nl ~ V. This implies that 
the lumpy structure should typically be saturated in q^'^-' 
with k < V/2. Consequently, g''^-' is expected to be well 
sufficient for all ensembles in Table 1. 

We now ask whether the subset UiCi of the lattice con- 
taining most of topological charge can be identified using 
q^*^^. A simple way to proceed is to order lattice points 
by the magnitude of qi'^'' and compute the running sum 
of positive {Q^''^{f)) and negative (Q~'*^(/)) charge as 
the fraction / of the highest points included increases. 
Functions Q^'*^ {Q~'^) should stabilize to a constant (« 
integer-valued) plateau at a well-defined value of / cor- 
responding to fraction of volume occupied by yJiCi. In 
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FIG. 2: The function Q+''=(/) for configuration C2 (Q = 1) of 
ensemble £5. The lowest curve corresponds to A; = 1. 




T average i — ■ — i 

1 _ \ config 1 

config 2 

I ' ' ' ' ' 

1 2 3 4 5 

V. 

FIG. 3: \Q\ +2k- Q'-'"*" for ensemble 64. 

Fig. Elwe show the behavior of Q^^'^ for configuration C2 
{Q = 1) from £5. We find no sign of plateaus for any fc, 
and nothing special happens around k — 2. Such behav- 
ior is characteristic for all configurations from the ensem- 
bles in Table n The same conclusion applies to functions 
Q'^if) monitoring total charge. The smooth monotonic 
behavior exhibited in Fig. |21 excludes the possibility that 
the bulk of topological charge is effectively concentrated 
in a small subvolume UiCi of typically isolated lumps. 

For another quantitative test, consider a configuration 
with topological charge Q and assume it is dominated by 
IQI -f ^ unit lumps (antilumps) and ^ antilumps (lumps). 
If the number of zeromodes is minimal, i.e. \Q\ = Nq 
(true for all our configurations), then the dimension of 
the topological subspace is = IQI + and we should 

have J2x k^^^l = Q^'"''" ~ Nl = \Q\ + 2C In fact, the 
quantity \Q\ + 2k — Q^^°-^'^ would be close to zero for 
all A: < ^ (since each mode in the topological subspace 
contributes approximately unity to Q^'"-^''), and it would 
start to increase rapidly for fc > ^. In a theory where 
the topological mixing scenario is relevant, monitoring 
the fc-dependence of IQj -I- 2A: — Q^'"-^^ could thus serve as 
a procedure for determining the dimension of the topo- 
logical subspace for a given configuration. We have com- 
puted this fc-dependence for all our configurations. The 
minimal value is always achieved at fc = 1 and is of order 
one rather than close to zero. Moreover, in ensembles 
£4, £5, where a wide range of fc is available, we observe 
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FIG. 4: Distribution of charges in coherent fluctuations for 
ensemble £a- 

monotonic (linear) increase for every configuration. This 
robust behavior is illustrated in Fig. O confirming again 
that there are no signs of a distinctive topological sub- 
space and hence no signs of dominance by unit lumps. 

3. The above arguments do not depend on specific 
properties of £i, such as their shape, volume or particu- 
lar field content. The results show that the topological 
charge is effectively carried by the bulk of the lattice, 
which is inconsistent with dominance by unit quantized 
lumps. At the same time, we find inhomogeneous behav- 
ior with noticeable peaks in g^*^-* , accompanied by visible 
coherence (see C^i^k) in Fig. ^| . As a first step toward un- 
derstanding this structure we now provide a simple char- 
acteristic of the typical values of charge associated with 
such coherent fluctuations. Given an arbitrary density q 
we consider the sets X' of centers of coherent behavior 
with lattice resolution namely 

= {x : \qy\ < \q^l q^qy >0 Wy; <\x^ y\^ <l} . 

The elements of X'' are local maxima oi \q\ over distance 
Vl, for which the sign of q is coherent over at least the 
same distance. Obviously, 1^ D 1^ D 1^ D . . .. To ev- 
ery X G X'' we assign a radius defined as the maximal 
distance from x over which the density is still coherent, 
i.e. Rx > Vl. The charge Qx is also assigned by sum- 
ming the density over the sphere of radius Rx centered 
at X. However, to be able to interpret Qx as a charge 
corresponding to an individual fluctuation, there should 
typically be no other centers within Rx. This allows for 
fixing the resolution in a self-consistent way by choosing 
the smallest I with this condition satisfied. 

To carry out the above procedure meaningfully for q^'^^ , 
one should work in the range of k where coherent fluctu- 



ations and their charges are relatively stable with respect 
to a change of k. For our ensembles, this is best satisfied 
for £4 where q^^-' is available. In physical terms, this cor- 
responds to including eigenmodes with imaginary parts 
up to « 500 McV. The fraction of overlapping centers 
drops dramatically (to 5 %) at / = 3, which is the value 
we have used. The corresponding distribution of charges 
for non-overlapping centers is shown in Fig. 0] A prop- 
erty which is insensitive to the choice of k and I is that 
the distribution effectively ends at about 0.5. This ap- 
pears to hold also in the case of full density qx (with only 
two configurations available) . Interestingly, this behavior 
might be compatible with the presence of center vortices 
in the QCD vacuum as the recent discussion of topology 
in the field of an idealized vortex suggests Ts] (The possi- 
ble manifestation of center vortices in topological charge 
fluctuations was also recently discussed in Ref. [T3|.'). 

To summarize, we have proposed that the low-energy 
behavior of topological charge density can be studied by 
using a suitable non-local realization of q{x) as a filter for 
short-distance fluctuations. If one is interested in the as- 
pects of topological charge affecting the low-momentum 
propagation of light quarks, then the appropriate non- 
local realization is naturally available through the low- 
eigenmode expansion of q{x). Such fermion filtering 
can be explicitly realized starting from lattice-regularized 
theory. We have proposed the expansion of qx associated 
with GW fermions as an ideal tool for this purpose. This 
provided us with a consistent framework to test whether 
the propagation of light fermion is effectively driven by 
the dominance of unit topological lumps in the QCD vac- 
uum. We find that this is not the case, implying that the 
picture of SxSB based on the mixing of corresponding 
topological "would-be" zeromodes is not accurate. We 
emphasize that we have not ruled out the logical possibil- 
ity that some unit quantized structures with well-defined 
boundaries can occur. However, we have shown that the 
hulk of topological charge does not come in that form. A 
first attempt to characterize the inhomogeneous nature 
of topological charge fluctuations in filtered densities re- 
sulted in very interesting results indicating that patterns 
of local behavior can provide us with detailed informa- 
tion about dynamically important structures in the QCD 
vacuum. 
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